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ABSTRACT 

Detailed  measurements  of  wall  pressure  fluctuations 
have  been  made  in  the  intermittent  region  of  a  flat  plate 
boundary  layer.  Digital  sampling  and  processing  techniques 
were  used.  The  properties  of  these  pressure  fluctuations  were 
found  to  be  similar  to  the  previous  measurements  made  in  the 
fully  turbulent  region.  The  measurements  were  repeated  with 
a  single  two  dimensional  surface  roughness  on  the  plate.  The 
only  changes  in  the  results  were  a  decrease  in  the  transition 
Reynolds  number  from  2  x  106  to  1.2  x  106  and  an  increase  in 
the  decay  rate  of  the  longitudinal  cross-spectral  density 
magnitude  by  a  factor  of  about  1.5.  Emmons'  analytical  model 
of  the  burst  rate  in  the  transition  region  was  found  to  be 
inaccurate.  His  model  treats  the  sources  of  the  turbulent 
spots  as  independent  random  events  with  prescribed  probability 
density  functions.  Both  a  delta  function  and  a  constant 
were  used  as  the  source  density  functions  and  in  each  case 
the  burst  rate  was  about  two  times  higher  than  the  present 
measurements . 
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1.  INTRODUCTION 

Wall  pressure  fluctuations  in  the  turbulent  boundary 
layer  have  long  been  recognized  as  important  sources  of  vib¬ 
ration  of  marine  and  aircraft  structures.  In  marine  applica¬ 
tions,  the  effect  of  turbulence  on  sonar  structures  is  of 
considerable  importance.  It  is  well  known  that  the  turbulent 
boundary  layer  is  a  major  source  of  sonar  self-noise. 

Boundary  layer  flows  have  received  much  attention  in  the 
past  and  an  excellent  list  of  publications  may  be  found  in 
White  [1],  Two  works  that  are  of  interest  in  this  study  are 
those  by  Blake  [2]  and  DeMetz  and  Casarella  [3],  Blake 
presented  a  thorough  experimental  study  of  the  wall  pressure 
fluctuations  in  fully  turbulent,  flat  plate  boundary  layers. 

The  statistics  he  measured  are  of  great  value  in  the  design  of 
structures  exposed  to  turbulent  flows.  DeMetz  and  Casarella 
investigated  the  intermittent  portion  of  a  flat  plate  boundary 
layer.  They  presented  many  details  of  the  intermittent  region 
but  did  not  measure  any  spatial  statistics.  The  properties  of  the 

intermittent  region  are  important  because  it  may  cover  most  of  a 

small  structure  and  cause  an  appreciable  amount  of  excitation. 

This  report  presents  the  results  of  an  experimental 
study  of  the  wall  pressure  fluctuations  in  the  intermittent 
region  of  transition  to  turbulent  flow  on  a  flat  plate. 

The  study  is  similar  to  the  one  of  DeMetz  and  Casarella  [3] 


but  includes  measurements  of  spatial  statistics  (cross- 
spectral  densities) .  Experiments  were  done  in  the  low  noise, 
low  turbulence  acoustic  wind  tunnel  in  the  Acoustics  and 
Vibration  Laboratory  at  the  Massachusetts  Institute  of 
Technology.  Two  phase-matched  condenser  microphones  with 
pinhole  caps  were  used  to  obtain  wall  pressure  signals.  The 
statistics  of  these  signals  were  calculated  digitally  using 
a  minicomputer.  The  intermittency ,  average  burst  rate  and 
average  burst  period  along  with  the  auto  and  cross-spectral 
densities  were  calculated  for  many  test  configurations. 
Convection  velocities  were  calculated  from  the  longitudinal 
cross-spectral  densities. 

The  results  obtained  compare  well  with  the  measure¬ 
ments  of  Blake  t2]  and  DeMetz  and  Casarella  [ 3 1 .  Most  of 
the  data  showed  no  dependence  on  intermittency.  The  auto- 
spectral  densities  varied  with  intermittency  but  this 
variation  was  most  likely  due  to  an  improper  scaling  along 
with  the  effect  of  spatial  averaging  over  the  microphone. 

The  measurements  were  repeated  with  a  single  down¬ 
stream  facing  step  in  the  plate  located  just  behind  the 
leading  edge.  The  results  were  very  much  the  same  as  for 
the  smooth  case.  The  only  noticeable  changes  were  an 
increase  in  the  decay  rate  of  the  longitudinal  cross-spectral 
density  magnitude  and  a  decrease  in  the  transition  Reynolds 


number. 
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2 •  PROPERTIES  OF  THE  TRANSITION  REGION 

In  the  flow  past  a  structure,  the  laminar  portion  of 
the  boundary  layer  is  usually  considered  passive  in  the 
excitation  of  the  structure.  This  includes  the  non-turbulent 
portions  of  the  transition  region.  However,  turbulent  spots 
can  contribute  an  appreciable  amount  of  excitation  and  must 
be  considered.  It  is  therefore  important  to  be  able  to 
predict  the  start  and  extent  of  the  transition  region  as 
well  as  its  important  statistical  properties.  The  transition 
process  has  received  much  attontion  in  the  past,  but  only  a 
small  number  of  investigations  have  been  done  on  the 
statistics  of  the  intermittent  region.  The  present  work  was 
concentrated  in  this  area. 

2 • 1  THE  TRANSITION  PROCESS 

The  natural  process  of  transition  to  turbulent  flow 
over  a  flat  plate  begins  with  the  instability  of  certain 
small,  random  velocity  fluctuations  in  the  laminar  flow. 
Fluctuations  with  certain  characteristics  become  unstable 
and  are  known  as  Tol lmien-Schlichting  (T-O)  waves.  These 
T-ft  waves  soon  develop  a  three  dimensionality.  This  is 
followed  by  vortex  breakdown  at  regions  of  high  localized 
shear  and  finally  the  formation  of  turbulent  spots  at  areas 
of  intense  fluctuations.  Turbulent  spots  (or  bursts)  are 
characterized  by  their  highly  random  velocity  fluctuations 
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and  small  scales  of  motion,  much  like  a  fully  turbulent 
boundary  layer.  Turbulent  spots  grow  as  they  are  convected 
downstream  and  soon  the  entire  boundary  layer  becomes 
turbulent.  Schubauer  and  Klebanoff  [4]  measured  the  size 
and  shape  of  a  typical  turbulent  spot  created  by  an  electric 
spark.  Their  findings  are  reproduced  in  Figure  2.1.  This 
figure  shows  the  normal  "arrowhead"  shaped  planform  growing 
at  a  half-angle  of  about  11°.  It  also  shows  a  profile  view 
of  the  spot.  Figure  2.2  is  an  idealized  sketch  of  the 
transition  process  taken  from  White  [l].  Transition  has 
received  much  attention  and  is  fairly  well  understood. 
Hundreds  of  studies  have  been  published  on  the  many 
interesting  and  often  complex  aspects  of  the  transition 
process . 

Emmons  [9]  suggested  that  burst  formation  be  modeled 
with  a  source  probability  distribution.  This  method 
requires  that  the  probability  of  a  turbulent  spot  forming 
at  a  certain  location  be  given  by  a  probability  distribution 
in  the  streamwise  coordinate.  Dhawan  and  Narasimha  [11] 
have  since  shown  experimentally,  that  burst  formation 
can  best  be  modeled  with  a  Dirac  delta 

function  probability  distribution.  The  location  of  the 
delta  function  is  usually  called  the  transition  point. 
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Dhawan  and  Narasimha  have  also  shown  that  the  boundary  layer 
velocity  profile  in  the  burst  region  (as  seen  from  a  fixed 
point)  fluctuates  between  a  Blasius  profile  between  bursts 
(originating  at  the  leading  edge)  and  a  turbulent  profile 
during  bursts  (beginning  at  the  transition  point).  Therefore, 
the  mean  velocity  is  a  time  average  of  the  laminar  and 
turbulent  portions.  This  observation  will  be  used  later  in 
Chapter  4. 

The  characteristics  of  wall  pressure  fluctuations 
are  of  interest  because  they  are  the  exciting  forces  in 
flow  related  structural  vibrations.  Typical  time  histories 
of  the  wall  pressure  in  the  transition  region  are  shown 
in  the  oscillographs  of  Figure  2.3.  A  turbulent  spot  appears 
as  a  high  frequency,  disordered  burst  in  the  otherwise  quiet, 
microphone  signal.  The  downstream  growth  and  convection 
of  the  spots  can  also  be  seen  in  this  figure.  The  lower 
signal  was  taken  at  a  position  slightly  downstream  of  that 
for  the  upper  signal.  A  turbulent  spot  in  the  upper  trace 
appears  slightly  larger  and  later  in  the  corresponding 
lower  trace  due  to  the  growth  and  convection. 

2.2  STATISTICAL  PROPERTIES 

The  one  most  important  statistical  property  of  che 
intermittent  region  is  the  intermittency  factor,  y  .  It 


Streamwise  Microphone  Separation 


Lateral  Microphone  Separation 


FIGURE  2.3  Time  Histories  of  Intermittent  Wall  Pressure 
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is  defined  as  the  fraction  of  time  the  flow  exhibits  a 
turbulent  nature  at  a  certain  point.  This  can  be  written: 

N 

Y  =  I  ti(/T,  (2.1) 

i=l  1 

where  t^  are  the  lengths  of  the  turbulent  time  periods,  and 
T  is  the  total  length  of  time.  Since  turbulent  spots  do 
not  have  square  profiles  (see  Figure  2.2),  y  varies  with 
vertical  as  well  as  horizontal  position.  Since  the  wall 
pressure  fluctuations  are  of  interest,  y  will  be  defined 
as  the  fraction  of  time  the  wall  pressure  exhibits  a 
turbulent  nature,  unless  it  is  specified  otherwise.  Notice 
that  this  is  not  the  same  as  y  defined  by  velocity 
fluctuations  near  the  wall  because  a  certain  amount  of 
averaging  takes  place  in  the  pressure  field  over  the 
microphone.  Obviously,  y  varies  from  0  (laminar)  at  the 
transition  point  to  1.0  (turbulent)  at  the  start  of  the 
fully  turbulent  boundary  layer.  It  should  be  pointed  out 
that  the  ends  of  the  transition  region  (y  =  0,  1.0)  are 
difficult  to  locate  accurately  due  to  the  random  nature  of 
the  transition  process.  Fair  approximations  can  be  made 
using  large  averaging  times. 

Two  other  statistical  properties  of  the  wall 
pressure  in  the  transition  region  are  the  average  burst 
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frequency  and  the  average  burst  period.  These  are  defined 
as  the  average  frequency  at  which  bursts  appear  and  the 
average  duration  of  the  bursts  appearing  at  a  point  in  the 
transition  region.  The  usual  units  are  burst  frequency  in 
bursts  per  second  and  burst  period  in  seconds.  It  can  be 
easily  demonstrated  that  of  these  three  statistical 
properties  -  intermittency,  burst  frequency  and  burst  period 
only  two  are  independent,  (i.e.  burst  period  can  be 
calculated  by  dividing  the  intermittency  by  the  burst 
frequency) . 

The  spectral  density  of  the  wall  pressure  is  an 
important  statistical  property.  It  is  very  useful  in 
determining  structural  excitation.  For  a  stationary, 
randomly  varying  wall  pressure,  p(t),  the  autocorrelation 
may  be  defined  by: 


Rp(t)  =  E[p(t)p(t  +  t)  ], 


(2.2) 


(see  Crandall  and  Mark  [ 12 J ) .  Here  E[  ]  denotes  an 
ensemble  average.  Taking  the  Fourier  transform  of  the 
autocorrelation  and  dividing  by  2tt  defines  the  spectral 
density  of  p(t),  4>p(u>): 


*p(w) 


-iwT  dx. 


Rp(t)  e 


(2.3) 
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Defined  in  this  manner  the  spectral  density  is  a  measure  of 
the  frequency  distribution  of  the  energy  in  p(t).  It 
can  be  shown  that  Rp  is  an  even  function  of  t.  Therefore 
^p(w)  is  an  even,  positive  real  function.  It  can  also  be  shown 
that  the  mean  square  pressure,  p2(t)  is  equal  to  the  area 
under  >^(10),  or: 

,  OO 

P2(t)  =  *  (uj)  dw.  (2.4) 

■'—CO  ^ 


The  spectral  density  is  a  good  measure  of  the  temporal 


characteristics  of  the  wall  pressure.  However,  it  provides 


no  information  regarding  the  spatial  variations.  One  method 
of  obtaining  spatial  information  is  to  make  two  wall  pressure 
measurements  simultaneously,  and  calculate  the  cross-spectral 
density.  This  method  was  used  by  Blake  [2],  Bull  [13], 
Willmarth  and  Wooldridge  [14 ]  and  others.  For  two  wall 
pressure  signals  separated  by  a  vector  distance,  r,  [p(x,t) 
and  p(x  +  r,  t) ]  the  cross-correlation  may  be  defined  by: 


Rp  ( r ,  t  )  =  E[p(x,t)p(x  +  r,  t+x)]. 


(2.3) 


E[  ]  again  denotes  an  ensemble  average,  and  p(x,t)  and 
p(x  +  r,  t)  must  be  stationary  in  time  and  spatially  in  the 
plane  of  the  wall.  The  cross-spectral  density  may  now  be 
defined  by: 
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^pCr.u))  *  |  Rp(r,T)e  iu'T  dT .  (2.6) 

Since  Rp(r,i)  is  not  necessarily  an  even  function,  4>(r,u>)  is 
in  general  complex.  It  is  usually  presented  as  a  magnitude 
and  a  phase.  The  magnitude  is  a  measure  of  the  correlation 
of  the  two  signals  at  each  frequency.  The  phase  can  serve 
as  a  measure  of  the  convection  velocity  of  turbulent  eddies. 
If  the  separation,  r  is  in  the  streamwise  direction,  the 
convection  velocity,  Uc  can  be  obtained  from  the  relation: 

Uc(tu)  =  uir/0  (u>)  (2. 7) 

where  0  (to)  is  the  phase  of  the  cross-spectral  density. 

Previous  studies  have  indicated  that  turbulent  spots 
have  many  of  the  same  characteristics  as  a  fully  turbulent 
boundary  layer.  Velocity  or  wall  pressure  signals  taken  in 
a  turbulent  spot  show  the  high  frequency,  apparently  random 
fluctuations  that  are  also  found  in  a  turbulent  boundary 
layer.  DeMetz  and  Casarella  [3]  have  shown  that  the  wall 
pressure  spectral  density  in  the  intermittent  region  has 
the  same  general  characteristics  as  a  turbulent  layer, 
although  the  magnitude  is  less  because  of  the  presence  of 
non-turbulent  regions  in  the  flow.  Cantwell,  et  al  [7] 
have  shown  that  the  boundary  layer  velocity  profile  in  a 
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turbulent  spot  contains  a  law  of  the  wall  renion 
and  a  wake  region,  just  as  is  found  in  a  turbulent  boundary 
layer.  These  observations  seem  to  indicate  that  a  turbulent 
spot  is  essentially  an  isolated  region  of  turbulence.  It 
would  be  helpful  if  some  of  the  statistics  in  the  inter¬ 
mittent  part  of  the  boundary  layer  could  be  compared 
directly  to  the  turbulent  data.  This  may  be  done  if  the 
intermittent  data  is  first  corrected  for  the  non-turbulent 
portions  of  the  flow.  DeMetz  and  Casarella  [3]  have 
shown  that  the  following  relationship  holds  for  the  mean 

square  wall  pressure  in  a  fully  turbulent  boundary  layer, 

~~2 

p  (t) i  and  the  mean  square  wall  pressure  in  an  inter- 

3 

mittent  boundary  layer,  p  (t)j^:* 

p2  (t)  |y  =  i  =  ~  P2(t)  |r  (2.8) 

Combining  equations  (2.4)  and  (2.S)  yields  the  relation: 

*p(u,,|Y=1  "  y  V")  lY*  (2-") 


♦Notice  that  p^TtT  | y  is  the  mean  square  of  the  total  pressure 
signal,  not  just  the  turbulent  part. 


1 
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Equation  (2.9)  allows  the  direct  comparison  of  the  turbulent 
statistic,  4>p  (u>)  |y  =  1  and  the  intermittent  one,  ♦pfwHy. 
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3.  EXPERIMENTAL  EQUIPMENT  AND  DATA  PROCESSING 

3.1  EXPERIMENTAL  FACILITY 

The  experiments  were  conducted  in  the  MIT  low- noise, 
low- turbulence  open  circuit  wind  tunnel.  This  facility  is 
shown  in  Figure  3.1  and  is  described  in  detail  in 
Hanson  [15].  It  has  a  flow  speed  range  of  20-50  m/s, 
controlled  by  varying  the  speed  of  a  DC  motor  driving  the 
blower.  Since  the  work  of  Hanson,  the  tunnel  has  been 
modified  to  improve  its  inherent  freestream  turbulence 
and  to  further  reduce  noise.  The  settling  chamber  has 
been  extended  and  the  number  of  turbulence  reducing  screens 
has  been  doubled.  Behind  the  screens  holding  the  honey 
comb  flow  straighener,  there  are  six  16  x  18  mesh  screens, 
followed  by  two  48  mesh  screens  and  a  two  foot  settling 
zone.  Figure  3.2  shows  the  1/3  octave  band  levels  of  the 
free  stream  turbulence  component  u'  obtained  for  various 
tunnel  speeds.  Figure  3.3  shows  the  resulting  percentage 
turbulence  level  as  a  function  of  free  stream  velocity. 

The  turbulence  level  is  approximately  0.05%.  Tani  [16] 
and  Spangler  and  Wells  [17]  have  shown  that  reducing 
free  stream  turbulence  below  this  level  does  not  change 
the  transition  Reynolds  number.  The  acoustic  noise  created 
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FIGURE  3.2  Ere  queries  Soectra  of  Free  Stream  Turbulence.  Measured 
at  Test  Section  Duct  Centerline  In  1/3  Octave  Bands 
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by  the  tunnel  has  been  reduced  by  the  isolation  of  the 
tunnel  and  blower,  the  installation  of  a  muffler-diffuser 
and  the  application  of  vibration  damping  and  sound  absorptive 
material  to  the  blower.  The  sound  pressure  spectrum  levels 
measured  just  outside  the  mixing  zone  in  an  open  jet 
configuration,  are  shown  in  Figure  3.4,  for  a  free  stream 
velocity  of  50  m/s.  It  is  felt  that  the  noise  levels  are 
low  enough, by  comparison  to  Spangler  and  Wells, as  to  not 
have  an  appreciable  effect  on  transition.  An  open  jet 
configuration  was  used  in  order  to  allow  some  of  the  sound 
generated  to  propagate  into  the  blockhouse.  The  block¬ 
house  was  airtight  and  completely  lined  with  two  inches 
of  polyurethane  foam  for  sound  absorption.  This  facility 
has  been  used  for  boundary  layer  research  in  the  past  by 
Blake  [2],  Martin  [18],  Shapiro  [19]  and  Moeller  [20]. 

3.2  TEST  FIXTURE 

Experiments  were  performed  on  a  flat  plate  mounted 
in  the  wind  tunnel  test  section.  The  test  section  was 
located  immediately  behind  the  contraction  and  just  inside 
the  blockhouse.  It  was  15  inches  (38  cm)  square  and  6  feet 
(183  cm)  long  and  was  lined  with  smooth  sound  absorptive 
material.  It  had  movable  walls  for  controlling  the 


FIGURE  3.4  Acoustic  Noise  in  Tunnel  Blockho 
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pressure  gradient.  The  test  section  was  also  equipped  with 
a  mechanical  x-y  traverse.  This  traverse  was  accurate 
to  0.001  inch  (0.025  mm)  vertically  and  to  0.1  inch  (0.25  cm) 
horizontally.  It  was  used  to  position  a  hot  wire 
or  a  pressure  probe  during  the  measurement  of  the  mean  flow 
properties  on  the  plate. 

The  test  plate  was  15  inches  (38  cm)  wide  and  48  inches 
(122  cm)  long  and  was  mounted  horizontally  in  the  test 
section  5  inches  (12.7  cm)  above  the  bottom  (see  Figures  3.5 
and  3.6).  The  plate  was  located  below  center  to  minimize  the 
secondary  flow  effects  described  by  Shapiro  [19].  The  plate 
was  constructed  as  a  composite  sandwich  to  increase  the 
vibration  damping  (see  Figure  3.7).  The  upper  layer  was 
made  of  0.25  inch  (0.635  cm)  aluminum  tool  and  jig  plate 
to  provide  a  flat,  easy  to  polish, surface.  The  center  layer 
was  made  of  a  0.125  inch  (0.317  cm)  thick  sheet  of  visco¬ 
elastic  vibration-damping  material  (E-A-R  C-1002,  from  E-A-R 
Corporation,  Westwood,  Massachusetts) .  The  lower  layer 
was  made  of  0.125  inch  (0.317  cm)  thick  aluminum  plate. 

The  whole  sandwich  was  bonded  together  with  a  two-part 
polyurethane  adhesive,  chosen  for  its  high  peel  strength. 

This  construction  technique  was  developed  by  Shapiro, 
and  resulted  in  a  damping  of  11  to  40%  of  critical  damping 


FIGURE  3.6  Tes 
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( depending  on  temperature  and  frequency) .  The  leading  edge 

of  the  plate  was  solid  0.50  inch  (1.27  cm)  thick  tool 

and  jig  plate  machined  into  a  6:1  ellipse  (see  Figure  3.8). 

This  shape  was  chosen  by  Shapiro  because  of  its 
practical  shape  and  lack  of  a  strong  pressure  gradient 
which  could  cause  premature  transition.  The  pressure 
transducers  were  flush  mounted  in  two  rotatable  3  inch 
(7.62  cm)  diameter  disks  (similar  to  DeMetz  and  Casarella's 
[3]).  Rotation  of  these  disks  allowed  transducer 
separations  of  0.275  to  5  inches  (0.7  to  12.7  cm,  see 
Figure  3.5).  The  trailing  edge  was  made  of  a  thin  splitter 
plate,  much  like  the  one  used  by  Shapiro,  which 
extended  rearward  to  prevent  coherent  vortex  shedding  from 
creating  plate  vibrations  (see  Figure  3.8). 

The  plate  was  assembled  and  the  upper  surface  was 
polished  to  a  surface  roughness  of  about  5u  inches  (0 . 127  microns) 
RMS.  It  was  then  mounted  in  the  test  section  and  adjusted 
to  have  a  minimum  pressure  gradient.  The  variation  in 
pressure  obtained  over  the  length  of  the  plate  at  a  free 
stream  velocity  of  35.8  m/s  is  shown  in  Figure  3.9.  The 
maximum  variation  was  about  1%  of  the  total  head  (excluding 
the  leading  edge) . 
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A  series  of  flow  visualizations  were  done  to 
investigate  the  effect  of  the  transverse  turbulence 
contamination  of  the  plate.  Figure  3.10  shows  the  extent 
of  the  contamination  for  a  free  stream  velocity  of  35.6  m/s. 
During  pressure  measurements,  the  transducers  were  always 
kept  at  least  0.48  inch  (1.2  cm)  from  the  edge  of  the 
contaminated  region  to  prevent  false  signals. 

Some  of  the  experiments  were  done  in  a  "rough" 
wall  configuration.  The  roughness  condition  consisted  of 
a  single,  downstream  facing  step  in  the  plate  (see 
Figure  3.11)  0.017  inch  (0.043  cm)  high,  and  4.25  inches 
(10.8  cm)  downstream  of  the  leading  edge.  This  type  of 
roughness  occurs  often  in  structures  at  welds  and  other 
seams.  Tani  [16]  determined  that  for  transition  to  occur 
at  such  a  roughness,  the  following  expression  should  be 
true : 


k  =  826.  (3.1) 

CO 

where  k  is  the  height  of  the  step.  In  the  present  work 
k  =  0.022  inch  (0.056  cm)  for  immediate  transition. 

Experiments  were  done  at  k  =  0.017  inch 
to  investigate  the  characteristics  of  transition  in  the  case 
where  it  begins  downstream  of  the  step. 


INSTRUMENTATION 


3.  3 

The  size  of  the  transition  region  (including  the 
starting  point  and  the  length)  were  measured  using  a 
hot  wire  anemometer  probe  (Thermo- Systems  Inc.  probe  model 
1261-T1,5)  with  a  Disa  model  55005  battery  powered, 
constant  temperature  anemometer  and  a  Disa  model  555D15 
linearizer.  A  schematic  of  this  equipment  is  shown  in 
Figure  3.12.  Static  pressure  measurements  were  made  with 
a  static  pressure  tube  and  a  Betz  micromanometer.  The 
resolution  of  this  micromanometer  was  0.1  mm  H20.  The  free 
stream  velocity  was  monitored  by  measuring  the  pressure  drop 
across  the  wind  tunnel  contraction  section  with  the 
micromanometer.  This  provided  a  resolution  of  about  0.1  m/s. 
A  calibration  was  performed  by  Shapiro  [19]  to  ensure  an 
accuracy  of  about  0 . 1  m/s . 

Wall  pressure  fluctuations  were  measured  using  two 
Bruel  andKjaertype  4138,  1/8  inch,  phase  matched  condenser 
microphones  fitted  with  1/32  inch  pinhole  caps.  The 
electronic  circuitry  is  shown  in  Figure  3.13.  Final 
statistical  calculations  were  made  digitally  using  an 
Interdata  model  M70  mini-computer  in  conjunction  with  a 
two  channel  eight-bit  analog  to  digital  converter.  The 
digital  calculations  were  made  using  the  data  acquisition 
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computer  program  described  below.  This  system  provided  a 
frequency  range  of  up  to  25,000  Hz  with  a  resolution  of 
1/256  of  the  maximum  frequency.  Phase  matched  amplifiers 
and  band  pass  filters  were  used  to  condition  the  signals 
prior  to  analysis  to  gain  the  maximum  signal  to  noise 
ratio  and  prevent  aliasing  problems.  A  single  channel 
spectrum  analyzer  was  used  to  detect  signal  contamination. 

By  placing  tape  over  the  pinhole  microphones  to 
eliminate  the  turbulent  pressure  fluctuations,  noise 
from  the  equipment,  from  vortex  shedding,  and  from  vibrations 
could  be  detected  and  dealt  with.  Figure  3.14  shows  tvpical 
spectra  of  the  taped  and  untaped  microphone  signals.  This 
figure  indicates  that  there  is  a  substantial  amount  of  noise 
below  1000  hertz.  This  noise  is  mainly  microphone  noise,  and 
acoustic  noise.  Measurements  were  made  with  an  accelerometer, 
which  indicated  extremely  low  vibration  levels.  Also  the 
coherence  between  a  pinhole  microphone  and  an  external 
microphone  was  measured.  These  signals  were  not  correlated, 
indicating  that  vortex  shedding  was  not  a  major  problem. 
Frequencies  below  1000  Hz  were  not  analyzed  to  avoid  noise 
contamination . 

The  calibration  equipment  is  also  shown  in  Figure  3.13. 
A  Oeneral  Radio,  model 1390-A  random  noise  generator  was  used  to 
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FIGURE  3.14  Turbulent  Pressure  Spectra  and  Noise  Floors 
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provide  a  wide  band  calibration  of  the  equipment  chain.  The 
results  of  this  test  are  shown  in  Figure  3.15.  The  equip¬ 
ment  sensitivity  was  determined  using  a  Bruel  and  Kjaer 
type  4220  pistonphone.  This  device  provided  a  sinusoidal 
pressure  fluctuation  of  known  magnitude  at  a  frequency  of 
250  Hz.  An  RMS  voltage  measurement  was  made  with  a  Bruel 
and  Kjaer  type  2607  measuring  amplifier.  The  sensitivities 
were  determined  as  channel  A:  -63.5dB  re  lV/pa, 
channel  B:  -65.4dB  re  lV/pa. 

The  microphones  were  tested  for  relative  phase 
match  by  Bruel  and  Kjaer  Instruments  Inc.  on  December  30,  1974. 
The  results  are  shown  in  Figure  3.16.  The  frequency  response 
of  each  microphone  was  also  checked,  using  a  Bruel  and 
Kjaer  type  4133  1/2  inch  condenser  microphone  as  a 
reference.  The  equipment  schematic  and  results  for  this 
test  are  given  in  Figures  3.17  and  3.18.  A  Helmholtz 
resonance  of  the  pinhole  cavity  was  noticed  at  about  17,000  Hz. 
Because  of  this  resonance,  the  analysis  was  extended  only  up 
to  13500  Hz. 

3.4  DATA  ACQUISITION  AND  PROCESSING 

The  data  acquisition  computer  program  was  divided 
into  three  major  parts,  (a  listing  can  be  found  in  Appendix  A). 
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FIGURE  3.15  Results  of  Equipment  Calibration 
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FIGURE  3.17  Microphone  Calibration  Set-up 
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FIGURE  3.18  Frequency  Response  of  Microphones 
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The  first  part  was  a  data  sampling  and  an  analog  to  digital 
conversion.  Then  a  conditioning  was  done  on  the  data  to 
determine  its  intermittent  characteristics.  The  method 
used  was  similar  to  that  of  Hedley  and  Keffer  [21], 

Finally,  a  spectral  analysis  was  done  on  the  conditioned 
signal. 

The  conversion  was  done  via  two  analog  to  digital 
(A/D)  converters  with  the  subroutine  CONVRT  (written 
by  G.  Holmes,  see  Appendix  A).  CONVRT  read  the  A/D 
converters  at  a  rate  equal  to  twice  the  maximum  frequency 
of  interest  to  prevent  aliasing.  This  data  was  stored  in 
a  data  array  for  further  processing. 

A  conditioning  was  then  done  on  the  entire  data 
array.  The  signal  was  first  squared  and  then  differentiated 
(using  a  central  difference  scheme) .  This  magnified  the 
inherent  differences  between  the  laminar  and  turbulent 
portions  of  the  signal.  The  signal  was  then  averaged  over 
successive  small  "smoothing  intervals"  to  smooth  out  sharp  fluc¬ 
tuations.  Finally,  the  conditioned  signal  was  compared  to 
a  predetermined  "criterion  level".  The  portions  of  time 
that  the  conditioned  signal  was  greater  than  the  criterion 
level  were  taken  as  turbulent  spots.  Since  the  character¬ 
istics  of  the  turbulent  spots  were  of  interest,  the  laminar 
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portions  of  the  signal  were  set  equal  to  zero.  This  reduced 
the  amount  of  non-turbulent  signal  present  (i.e.  Tollmien- 
Schlichting  waves,  noise,  etc.).  This  entire  conditioning 
process  is  shown  schematically  in  Figure  3.19. 

The  conditioning  process  had  two  variables,  smoothing 
interval  and  criterion  level.  These  variables  were  adjusted 
so  that  small  changes  in  them  would  not  cause  noticeable 
changes  in  the  resulting  conditioned  signal.  Hedley  and 
Keffer  [21]  have  shown  that  the  burst  rate  is  the  statistic 
most  sensitive  to  changes  in  the  conditioning  of  the  signal. 
Figures  3.20  and  3.21  show  the  effect  of  smoothing  interval 
and  criterion  level  on  the  burst  rate.  Figure  3.22  shows 
a  typical  raw  signal  with  a  "criterion  function"  superimposed 
on  it.  The  criterion  function  was  equal  to  one  during 
turbulent  spots  and  zero  during  the  laminar  portions  of  the 
signal.  This  figure  provides  a  good  check  of  the  signal 
conditioning. 

The  final  part  of  the  data  acquisition  computer 
program  was  to  perform  the  statistical  calculations.  First 
the  intermittency  function,  burst  rate  and  burst  period  were 
calculated  for  the  data  array.  Then  a  spectral  analysis 
similar  to  that  of  DeJong  [22]  was  performed  on  the  data 
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FIGURE  3.20  Burst  Rate  Versus  Smoothing  Interval 
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FIGURE  3.21  Burst  Rate  Versus  Triqqering  Level  (Criterion 
Level  ) 
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The  spectral  analysis  involved  calculating  spectral 
estimates  from  the  random  pressure  signals.  The  theory 
of  random  signal  processing  is  well  developed  (see  Bendat 
andPiersol  [23]  or  Oppenheim  and  Schafer  [24]).  The  wall 
pressure  spectral  density  was  defined  in  Chapter  2  by: 

i  f00 

V00  =  27  j  vT>e~la,TdT'  (2,3) 

where  Rp(r)  was  defined  as: 

Rp(t)  =  E[p(t)p(t  +  t)],  (2.2) 

and  p(t)  must  be  stationary.  If  p(t)  is  also  ergodic,  (2.2) 
may  be  replaced  by: 

(t)  =  lim| 

T-H» 

=  <p(t)p(t  +  t) >  (3.2) 


T/2 

p(t)p(t  +  T )  dt 

-T/2 


(<  >  denotes  a  temporal  average) .  The  assumption  of  station- 
arity  is  good  if  the  statistics  of  the  pressure  signals  do 
not  change  over  relatively  long  periods  of  time.  Although 
it  cannot  be  rigorously  proven,  the  signals  were  assumed  to 
be  ergodic.  This  allowed  the  temporal  average  to  be  used, 
which  provided  an  easier  calculation  method.  The  Fourier 
transform  is  used  extensively  in  the  analysis  of  random 
signals.  A  useful  approximation  is  the  finite  Fourier 
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P  is  a  function  of  the  integration  interval,  T  and  the 
frequency,  f  (f  =  u)/2v) .  p(t)  must  be  stationary.  If  the 
assumption  of  ergodicity  is  correct,  the  finite  Fourier 
transform  and  the  temporal  average  may  be  used  to  calculate 
the  spectral  density  instead  of  equation  (2.4). 

♦_(f)  =  lim  ~  < | P (f ,T) | 2>  (3.4) 

‘  T~*oo 

Notice  that  since  the  experimental  frequency,  f,  is 
being  used, 

V«>  '  3?  Vf)  <3-5> 

The  finite  Fourier  transform  may  be  approximated  by  the 
sum: 

N-l  .  . 

P(f,T)  =  h  £  p (nh) e”1<J7Trn  ,  (3.6) 

n=o 

if  the  function  p(t)  is  sampled  N  times  with  a  time  interval 

h.  Here  T  is  equal  to  Nh.  An  efficient  method  of  computing 

this  sum,  is  the  fast  Fourier  transform,  or  FFT ,  developed 

by  Cooley  and  Tukey  [25].  The  FFT  requires  only  2N  log2  N 

2 

multiplications,  instead  of  N  ,  and  provides  the  value  of 
P(f,T)  at  frequencies  of  fj  *  j/T.  The  FFT  requires  N  to 
be  an  integral  power  of  2  and  the  result  may  be  written  as: 
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Pj  =  (f  j  ,T)  ,  j=0,l,...|  -  1. 

An  estimate  of  4>  (f)  may  now  be  written  as: 

r 

VV  =  IT  ^jl2'  3  =  0,1,. ..f  -  l. 


(3.7) 


(3.8) 


Due  to  the  randomness  of  the  function  p(t),  the  standard 
error  of  the  estimate,  *p(f^)  is  1.  If  ip(fj)  is  calculated 
for  a  large  number  of  data  sets,  the  standard  error  may  be 
reduced  to  /l/M  (where  M  is  the  number  of  data  sets)  by 
taking  the  average,  ♦p(fj),of  all  the  4>p(fj)'s.  We  now 
have : 


VV  -ft  £  Vk'V 


F  (f.)  =  i  "  2H 
P*V  M  N 


M  N  lPjl  '  ^  ~  0,1, ...j  -  1.  (3.9) 


If  two  random  pressure  signals,  p(t)  and  q(t),  are 
sampled  simultaneously,  a  joint  spectral  analysis  may  be 
done  using  the  theory  of  complex  numbers.  If  a  complex 
number  is  defined  as: 


z(t)  =  p(t)  +  i  q(t)  , 


(3.10) 


Then  a  complex  Fourier  transform  may  be  approximated  by: 


Z.  ■  E  [p(nh)  +  i  q  (nh)  ]< 


i2Ijn 

N 


j  =  0,1, ...f  -  1 


(3.11) 


This  relation  now  provides: 


?  Czi  *  K-i> 


■k 


(3.12) 


where  *  denotes  complex  conjugate.  4>p(f^)  and  ^(f^)  way 
now  be  obtained  from  these  relations  and  equation  (3.8). 
The  cross-spectral  density  is  defined  as: 


Jr,f)  =  lim  E4rP(f,T)Q*(f,T)]. 


(3.13) 


r  is  tne  vector  between  pressures  p(t)  and  q(t). 
4>p(r,f)  may  be  approximated  by: 


y?'V  Sh  pj  °3*' 


(3.14) 


The  standard  error  of  4^r ,  )  can  be  reduced  by  averaging, 

as  with  $p. 

The  data  acquisition  computer  program  performed 
these  complex  FFT  calculations,  obtaining  the  pressure 
spectral  densities  and  the  cross-spectral  density.  All 
of  the  spectra  were  multiplied  by  1/y  to  correct  for  the 
intermittency.  During  the  data  acquisition  M  was  set  at 
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250  so  that  the  standard  error  was  6.3%. 

Another  computer  program  was  used  to  non-dimensiona- 
lize  and  plot  the  spectral  data  and  to  calculate  the 
convection  velocity.  Both  of  these  programs  may  be  found 
in  Appendix  A. 
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4.  DISCUSSION  OF  RESULTS 

4.1  MEAN  VELOCITY  MEASUREMENTS 

The  characteristics  of  the  mean  velocity  profiles  in 
the  intermittent  region  are  listed  in  Table  4.1  for  various 
flow  speeds.  In  the  intermittent  region  the  boundary  layer 
had  two  distinct  behaviors  -  laminar  and  turbulent.  Dhawan 
and  Narasimha  [ll]  have  shown  that  mean  velocity  measurements 
in  this  region  result  from  time  average  over  a  laminar  profile 
and  over  a  turbulent  profile.  They  have  also  shown  that  a 
very  good  approximation  of  the  flow  characteristics  are 
obtained  by  averaging  a  Blasius  boundary  layer  with  a  1/7  power 
law  turbulent  boundary  layer.  However,  a  virtual  origin  at 
the  transition  point  must  be  used  for  the  turbulent  boundary 
layer.  Figure  4.1  shows  a  typical  mean  velocity  profile  that 
we  measured  in  the  intermittent  region.  A  law  of  the  wall 
profile  and  a  Blasius  profile  are  also  shown  for  comparison. 

The  measured  displacement  thickness,  6*,  can  be 
expressed  by: 

6*  -  Y6J  ♦  (1  -  y)S*  ,  (4.1) 

where  5£  is  the  displacement  thickness  of  the  Blasius  boundary 
layer  and  is  that  for  the  turbulent  layer.  A  similar 
expression  can  be  written  for  the  coefficient  of  skin  friction, 

Cfm!  Cfm  =)Ca  +  (l-Y)CfT.  (4.2) 

Since  it  is  the  behavior  of  the  flow  in  the  turbulent  spots 
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TABLE  4.1 


6*Txl03 

cfTxio3 
(x=0. 91m) 

Velocity 

X 

-< 

il 

o 

• 

o 

Xl-y=0.99 

(x=0. 91m) 

(m/s) 

(m) 

(m) 

Eq .  4.7 
(m) 

Eq.  4.3 

Smooth  Condition 

36.0 

0.41 

1.09 

1.22 

3.65 

37.6 

0.36 

1.02 

1.32 

3.56 

38.3 

0.33 

0.94 

1.37 

3.51 

39.0 

0.30 

0.89 

1.43 

3.47 

39.5 

0.28 

0.81 

1.47 

3.44 

40.0 

0.25 

0.76 

1.53 

3.41 

40.5 

0.23 

0.71 

1.  56 

3.38 

44.5 

0.20 

0.66 

1.60 

3.30 

45.0 

0.18 

0.63 

1.64 

3.28 

45.  5 

0.17 

0.62 

1.65 

3.27 

Rough  Condition 

22.0 

0.64 

1.07 

0.77 

4.49 

22.5 

0.61 

1.02 

0.84 

4.38 

23.0 

0.57 

0.99 

0.94 

4.26 

24.0 

0.39 

0.94 

1.34 

3.90 

25.0 

0.29 

0.79 

1.55 

3.75 

25.5 

0.25 

0.71 

1.63 

3.69 

26.0 

0.23 

0.64 

1.66 

3.66 

26.5 

0.20 

0.53 

1.72 

3.62 

29.0 

0.17 

0.33 

1.76 

3.  53 

29.5 

0.15 

0.30 

1.80 

3.51 
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that  xs  o t  interest,  S*  and  CfT  are  the  desired  parameters. 

Measurements  of  skin  friction  were  not  made,  but  CfT  was 

determined  from  the  approximate  formula : 

CfT  =  0.455/fcn2  (0.06  Re(x-xt»,  (4.3) 

(White  Cl])-  The  Reynolds  number  is  based  on  (x-xfc)  since 

a  virtual  origin  at  the  transition  point  (xfc)  was  assumed. 

The  actual  transition  point  was  not  located  precisely,  but 

it  was  approximately  determined  as  the  point  at  which  the 

intermittency  was  0.01.  The  hot  wire  anemometer  and  an 

oscilloscope  were  used  to  locate  this  point. 

6*  and  9  were  determined  using  trapezoidal  rule 
mm 

integrations  of  the  integrals  in  their  definitions. 


(4.4) 


(4.5) 


where  <$  is  the  boundary  layer  thickness.  6*  was  determined 
using  the  following  formula  from  Blasius'  theory: 


6*  =  1.7208x//Re“  .  (4.6) 

li  X 

6*  was  then  calculated  using  equation  4.1.  For  this  cal¬ 
culation,  y  was  measured  using  the  data  acquisition 
computer  program  with  the  hot  wire  probe  again  placed 
a  distance  of  about  46*  above  the  plate.  6,f  was  also  calculated 
using  the  approximation: 
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(4.7) 


(White  [1]).  The  measured  values  of  6*  are  plotted  in 
Figure  4.2  against  those  calculated  from  equation  4.7. 

4.2  MEAN  PRESSURE  MEASUREMENTS 

Table  4.2  gives  a  summary  of  the  mean  and  intermittent 
characteristics  of  the  wall  pressure  field.  The  extent  of 
transition  on  the  test  plate  is  shown  in  Figure  4.3.  The 
variation  of  intermittency  with  Reynolds  number  (based  on 
x,  the  distance  from  the  leading  edge  of  the  plate)  is  shown 
for  both  the  smooth  and  rough  conditions.  Recall  that  in 
the  rough  condition,  the  test  plate  had  a  0.015  inch 
(0.381  mm)  high  two  dimensional  reverse  step  located  4.25 
inches  (10.8  cm)  behind  the  leading  edge. 

In  both  cases,  the  variation  of  intermittency  is  in 
the  form  of  a  Gaussian  integral  curve.  Similar  results 
have  been  found  experimentally  by  Schubauer  and  Klebanoff  [4] 
and  DeMetz  and  Casarella  [3]  and  analytically  by  Emmons  [9]. 
The  effect  of  the  roughness  was  to  reduce  the  transition 
Reynolds  number  from  about  2  x  106  to  1.2  x  10^. 

The  non-dimensional  burst  frequency  (fB*)  is  plotted 
against  intermittency  in  Figures  4.4  and  4.5.  Theoretical 
results  and  the  experimental  result  of  Farabee ,  et  al  [10] 
are  also  shown.  The  theoretical  results  were  obtained  by 
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TABLE  4.2 


Run  No. 

Velocity 

(m/s) 

X 

(m) 

Y 

.Bursts. 

sec 

tb 

(msec) 

Smooth 

Condition 

3606 

36.0 

0.91 

0.062 

73.2 

0.85 

3824 

38.3 

0.91 

0.245 

60.8 

4.03 

3825 

37.6 

0.91 

0.254 

214.0 

1.19 

3940 

39.0 

0.91 

0.399 

73.7 

5.42 

3864 

38.3 

0.91 

0.  644 

69.  9 

9.22 

4069 

39.5 

0.91 

0.688 

66.3 

10.4 

4089 

39.5 

0.91 

0.894 

32.2 

27.8 

4076 

40.5 

0.91 

0.765 

64.4 

11.9 

4499 

44.5 

0.91 

1.000 

0.316 

3160 

4599A 

45.0 

0.91 

0.997 

1.  50 

665 

3818 

37.6 

0.91 

0.184 

48.7 

3.78 

3922 

39.0 

0.90 

0.218 

52.3 

4.16 

3840 

38.3 

0.91 

0.397 

71.1 

5.58 

4041 

39.5 

0.90 

0.406 

70.8 

5.74 

3954 

39.0 

0.91 

0.539 

77.4 

6.96 

4046 

40.0 

0.90 

0.464 

71.9 

6 . 46 

4078 

40.0 

0.91 

0.775 

59.2 

13.1 

4065 

40.5 

0.90 

0.652 

76.7 

8.50 

4599B 

45.  0 

0.91 

0.998 

1. 18 

843 

4699 

45.5 

0.90 

0.996 

2.53 

394 

Rough 

Condition 

2219 

22.0 

0.91 

0.188 

43.7 

4. 31 

2322 

23.0 

0.91 

0.219 

46.7 

4.70 

2326 

23.0 

0.91 

0.260 

56.0 

4.65 

2435 

24.0 

0.91 

0.354 

70.8 

4.94 

2440 

24.0 

0.91 

0.401 

76.7 

5.23 

2441 

24.5 

0.91 

0.405 

86.1 

4.71 

2570 

25.0 

0.91 

0.701 

81.4 

8.62 

2670 

25.5 

0.91 

0.698 

72.1 

9.67 

2999A 

29.0 

0.91 

0.996 

5. 13 

197 

3099 

29.5 

0.91 

0.993 

7.90 

126 

2225 

22.5 

0.91 

0.247 

55.4 

4.46 

2325 

23.0 

0.90 

0.251 

51,5 

4.88 

2448 

24.0 

0.91 

0.475 

80.4 

5.91 

2639 

25.5 

0.90 

0.389 

76.7 

5.07 

2562 

25.0 

0.91 

0.624 

84.9 

7.35 

2652 

26.0 

0.90 

0.525 

82.8 

6.34 

2685 

26.0 

0.91 

0.847 

53.7 

15.  8 

2674 

25.5 

0.90 

0.734 

70.9 

10.4 

2999B 

29.0 

0.91 

0.989 

9.24 

10.7 

2690 

26.5 

0.90 

0.905 

38.9 

23.  3 

applying  Emmons'  [9]  theory  to  two  different  source  rate 
density  functions.  In  Figure  4.4  a  line  source  density 
function  was  assumed  and  a  constant  source  density  function 
was  used  in  Figure  4.5.  In  the  constant  case,  the  source 
density  function  was  constant  downstream  of  the  transition 
point,  but  was  set  to  zero  upstream.  In  the  line  source 
case  Dirac's  delta  function  was  used  as  the  source  rate 
function.  The  delta  function  was  centered  on  the  transition 
point.  The  non-dimensionalization  was  done  as  in  Farabee,  et 
al.  For  the  line  source: 

fB*  E  fB/n  tana  '  (4*8) 

and  for  the  constant  source: 


fB/g  tana 


2/3 


(4.9) 


where  ffi  is  the  burst  frequency,  n  is  the  line  source  density 
and  g  is  the  uniform  source  density.  The  burst  growth  angle 
is  a  and  a*  is  a  non-dimensional  parameter  defined  by, 

U  Uoo 

a*  =  tana  .  (4.10) 

U*Ut 

where  is  the  velocity  of  the  leading  edge  of  the  burst, 

Ut  is  the  trailing  edge  velocity  and. 


-72- 


The  numerical  values  for  V z>  Ufc  and  a  were  taken  from 
Schubauer  and  Klebanoff  [4]  (see  Figure  2.1). 

If  a  line  source  function  is  used  in  Emmons'  model 
of  the  intermittent  region,  the  following  relations  are 
obtained. 


fB*  =  2  J ( 1-Y )  £n  (I~-) 


and 


V  =  0.420  fB 


AX„ 


U 


where 


AxT  Xly=.99  “ 


Y  =  .  01 


(4.12) 


(4.13) 


(4.14) 


The  use  of  a  constant  source  function  gives  the  relations: 


ffi*  =  8.722  [  nexp[3.035(n3-x3) ]dn, 

'  n 


(4.15) 


and 


where 


fB*  =  0.412 


AX_ 


•B 


U. 


(4.16) 


(4.17) 


Equations  4.12,  4.13,  4.15  and  4.16  are  obtained  in  Appendix  B. 
Equations  4.12  and  4.15  were  used  to  obtain  the  theoretical 
curves  in  Figures  4.4  and  4.5.  The  experimental  results  were 
scaled  using  equations  4.13  and  4.16. 
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The  burst  rate  measurements  of  Farabee,  et  al  [10] 
compare  well  with  the  present  results  for  the  smooth  plate. 

The  rough  plate  results  are  higher  than  the  smooth  case,  but 
this  is  most  likely  due  to  inaccuracies  in  the  measurement 
of  AXt.  The  burst  rates  were  measured  with  the  pin-hole 
microphones  and  AXT  was  measured  with  the  hot  wire  anemometer 
at  a  later  time.  If  the  source-rate  function  changed  during 
this  time,  the  measured  value  of  AxT  would  be  in  error. 

The  theoretical  curve  of  Figure  4.4  (line  source  model) 
has  approximately  the  same  shape  as  the  measured  results 
but  its  magnitude  is  too  high.  Emmons  [9]  points  out  that 
this  theory  will  predict  higher  burst  rates  than  are  measured. 
When  two  bursts  are  separated  by  a  very  small  laminar  region, 
a  transducer  may  only  detect  one  burst  whereas  the  theory 
will  count  two  separate  burst.  (This  is  due  to  the  smoothing 
in  the  data  acquisition  computer  program  . )  If  two  or  more 
bursts overlap,  the  theory  treats  them  as  a  single  burst  and 
uses  the  period  of  the  larger  burst.  However,  if  the  bursts  do 
not  completely  overlap,  the  observed  period  is  somewhat 
longer  than  this.  These  effects  tend  to  make  the  theoretical 
burst  rate  higher  than  the  observed  rate.  A  basic  assumption 
in  Emmons'  theory  is  that  the  turbulent  spots  are  completely 
independent  of  each  other.  This  assumption  is  not  totally 
correct  but  seems  to  be  a  good  approximation,  since  strong 
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burst  interaction  would  change  the  shape  of  the  curve  in 
Figure  4.4. 

The  mean  square  pressure  ratio  (p2 | y/p2 I y=\)  is 
plotted  against  intermittency  in  Figure  4.6.  The  mean  square 
pressure  ratio  is  defined  as  the  ratio  of  the  total  mean 
square  pressure  to  the  mean  square  of  the  turbulent  portion 
only.  If  the  laminar  portions  of  the  intermittent  signal 
do  not  contribute  to  the  mean  square  pressure  then  the 
following  relation  from  Section  2.2  must  be  true: 


The  mean  square  pressure  ratio  must  then  be  equal  to  y.  Due 
to  noise  and  laminar  pressure  fluctuations  (e.g.  T-S  waves) 
the  mean  square  pressure  ratio  was  higher  than  y  when  y 
was  less  than  about  0.6.  The  difference  increased  for 
smaller  y  because  of  the  greater  amount  of  laminar  signal. 
4.3  STATISTICS  OF  PRESSURE  FLUCTUATIONS 

The  wall  pressure  spectral  densities  for  the  smooth 
and  rough  conditions  are  shown  in  Figures  4.7  and  4.8.  The 
spectral  density  was  non  dimensionalized  by  Uoo/q26*y,  where 
q  is  the  total  head  (1/2  pu* ) ,  and  the  frequency  was  scaled 
by  6^/U^.  Results  of  Blake  [2]  and  Willmarth  and  Wooldridge 
[14]  for  fully  turbulent  boundary  layers  are  shown  for 
comparison.  The  spectral  density  is  shown  for 
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several  values  of  intermittency  (y)  ,  ranging  from  0.06  to 
0.99.  Each  curve  has  been  faired  through  201  data  points. 

The  maximum  scatter  in  the  data  was  1  dB.  The  figures  show 
that  the  magnitude  of  the  wall  pressure  spectral  density 
increased  by  about  10  dB  as  y  varied  from  0.99  to  0.06.  The 
results  of  Blake  and  Willmarth  and  Wooldridge  are  about  3  dB 
below  the  y  =  0.99  case.  These  variations  in  magnitude  are 
most  likely  due  to  an  improper  scaling  of  the  spectral  density. 
In  the  frequency  range  of  the  present  measurements  (0.2  < 

<  5)  the  wall  shear  stress  of  the  turbulent  part  of 
the  flow  (t^)  is  a  better  parameter  than  q  because  it  is 
related  to  the  flow  behavior  near  the  wall,  where  the  pressure 


L 


sources  are  located  (see  Blake).  Figure  4.9  shows  the  wall 


pressure  spectra,  rescaled  by  U0o/Tw^,<5*y  along  with  the  results 


of  Blake.  A  much  better  collapse  of  the  results  is  achieved. 
The  turbulent  wall  shear  stress  was  calculated  using 

0.455 


wT 


-fT 


q  en2(0..06Re.  ,) 

(x-xt> 


(4.3) 


(White  [1]) . 

The  low  intermittency  results  are  not  shown  in 
Figure  4.9  because  equation  4.3  does  not  give  accurate  results 
if  x-xt  is  small.  This  equation  was  originally  obtained  for 
fully  turbulent  boundary  layers  and  is  not  accurate  at  low 
Reynolds  numbers.  Also,  the  error  in  (x-xfc)  was  large  when 
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x  was  near  xfc  due  to  the  inaccuracy  in  xfc.  The  transition 
point  (actually  x/^_Q  was  measured  after  the  pressure 
fluctuations  were  measured  so  that  slight  changes  in  the  flow 
characteristics  could  have  caused  appreciable  errors  in 

The  spectra  in  Figures  4.7  and  4.8  roll  off  faster 
at  the  low  intermittencies  than  at  the  higher  values.  This 
is  due  to  the  effect  of  spatial  averaging.  The  reduced 
frequency  (ukS*,/Uiv)  at  which  spatial  averaging  is  important 
depends  on  the  ratio  d/<S*,  where  d  is  the  microphone  pinhole 
diameter.  The  ratio  d/<S *i  is  given  for  each  curve  in  the 
figures.  The  lower  intermittency  cases  had  larger  values  of 
d/iS^,  so  that  their  spectra  roll  off  at  lower  reduced  frequencies 
than  for  the  higher  intermittencies.  In  Figure  4.8,  d/<5* 
varies  from  0.48  to  0.99  and  the  change  in  the  roll-off 
frequency  is  most  noticable. 

The  wall  pressure  convection  velocities  for  the 
smooth  and  rough  conditions  are  shown  in  Figures  4.10  and 
4.11.  Each  curve  in  these  figures  has  been  faired  through 
201  data  points,  all  taken  at  a  specified  intermittency.  The 
■4C.it ter  in  the  data  for  each  curve  was  about  10%.  The  fully 
* urbulent  results  from  Blake  [2]  and  Bull  L 1 3 ]  are  shown 
xnparison.  The  agreement  with  the  phase  convection 
!.•»  of  Blake  is  good.  This  is  to  be  expected  since  a 
•  >>  ‘w  ith  was  used,  whereas  Bull's  group  velocities 
.  i  with  a  1/3  octave  filter.  These  results 


Uc  /  U 


indicate  that  convection  velocity  does  not  depend  on  inter- 
mi  ttency  and  that  the  single  roughness  has  no  effect.  The 
y=0.19  curve  in  Figure  4.11  goes  only  to  ^6*/^  =  1.0.  This 
is  because  the  magnitude  of  the  wall  pressure  spectrum  dropped 
abruptly  above  this  frequency  and  the  convection  velocity 
could  not  be  calculated. 

Figures  4.12  and  4.13  show  the  normalized  longi¬ 
tudinal  cross-spectral  density  magnitudes  of  wall  pressure 
for  the  smooth  and  rough  conditions.  Again  the  curves  in 
these  figures  have  been  faired  through  201  data  points.  The 
scatter  in  the  data  was  about  10%.  The  cross-spectral  density 
magnitude  was  divided  by  the  spectral  density  and  the 
frequency  was  non-dimensionalized  by  rg/Uc,  where  rs  is  the 
streamwise  microphone  separation  and  Uc  is  the  phase  convection 
velocity.  These  curves  show  the  coherence  in  the  pre¬ 
ssure  field  versus  the  microphone  separation  in  wave  length 
units.  The  results  of  Blake  [2]  and  Willmarth  and  Wooldridge 
[14]  are  shown  for  comparison.  Our  results  do  not  indicate 
any  dependence  on  intermi ttency .  However,  the  curves  are 
slightly  lower  than  the  results  of  Blake  and  Willmarth  and 
Wooldridge.  This  difference  can  best  be  explained  by  the 
fact  that  in  the  intermittent  region,  a  turbulent  spot  will 
not  cover  both  microphones  all  of  the  time.  Thus,  for  part 
of  the  time  one  microphone  will  have  a  laminar  signal  while 
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FIGURE  4.12  Normalized  Lonqitudinal  Cross-Spectral  Densities  (Smooth  Wall). 
Y  is  shown  in  percent. 


FIGURE  4.13  Normalized  Longitudinal  Cross-Spectral  Densities  (Rough  Wall) 
Y  is  shown  in  percent. 
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the  other  has  a  turbulent  signal  (see  Figure  2.3).  This 
reduces  the  magnitude  of  the  cross-spectral  density.  The 
results  for  the  rough  wall  were  similar  to  those  for  the 
smooth  wall,  but  decayed  1.5  times  faster.  Blake's  rough 
wall  curve  decayed  much  faster  but  it  was  for  a  more  severe 
roughness  condition. 

The  normalized  lateral  cross-spectral  density 
magnitudes  of  the  wall  pressure  for  the  smooth  and  rough 
conditions  are  shown  in  Figures  4.14  and  4.15.  The  results 
for  the  fully  turbulent  boundary  layers  of  Blake  [2]  and 
Bull  [13]  are  also  shown.  The  vertical  bars  show  the  maximum 
variation  in  the  data  for  all  intermittencies  measured. 
Measurements  were  made  at  five  values  of  intermittency  for 
each  figure  ranging  from  0.28  to  0.99.  There  were  no 
variations  in  the  magnitude  which  could  be  attributed  to 
changes  in  intermittency.  These  results  compare  well  with 
the  results  of  Blake  and  Bull. 
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FIGURE  4.14 


Normalized  Lateral  Cross-Spectral  Densities  (Smooth  Wall). 
Y  is  shown  in  percent. 
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FIGURE  4.15  Normalized  Lateral  Cross-Spectral  Densities  (Rough  Wall) 
Y  is  shown  in  percent. 
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5.  CONCLUSIONS 

Measurements  of  the  statistics  of  the  wall  pressure 
fluctuations  in  the  intermittent  region  of  a  viscous  boundary 
layer  on  a  flat  plate  have  been  presented.  They  were  found 
to  be  similar  to  the  results  obtained  in  the  fully  turbulent 
region  by  Blake  [2],  Bull  [13],  Willmarth  and  Wooldridge  [14] 
and  others.  When  properly  scaled,  the  statistics  agreed  with 
the  results  cited  in  previous  reports  for  the  fully  turbulent 
region  (to  experimental  accuracy) .  The  wall  pressure  spectral 
densities  did  not  all  collapse  to  the  same  curve,  but  this 
was  due  to  the  inaccurate  determination  of  the  wall  shear 
stress  and  the  effect  of  spatial  averaging  over  the  microphone 
pinhole.  The  wall  shear  stress  is  essential  to  the  correct 
scaling  of  the  spectral  density.  The  reduced  frequency 
(coS^/U^)  at  which  spatial  averaging  becomes  important  depends 
on  the  ratio  of  the  microphone  pinhole  diameter  to  the  tur¬ 
bulent  boundary  layer  displacement  thickness.  The  similarities 
in  the  properties  of  turbulent  spots  and  turbulent  boundary 
layers  seem  to  indicate  that  the  underlying  physical  phenomenon 
is  the  same  in  each  case. 

The  structural  excitation  due  to  the  intermittent 
region  of  a  boundary  layer  may  therefore  be  determined  by 
using  results  for  fully  turbulent  layers  if  the  effect  of 
the  intermittency  is  accounted  for.  The  location  of  the 
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transition  point  and  the  variation  of  intermittency  must 
be  estimated  first.  These  calculations  should  include  the 
influence  of  surface  roughness,  free  stream  turbulence  and 
other  external  effects.  The  variations  in  intermittency 
measured  here  were  found  to  be  similar  to  previous  analytical 
and  experimental  results  (Emmons  [9],  DeMetz  and  Casarella  [3] 
and  Schubauer  and  Klebanof f  [4 ])  . 

The  burst  frequency  was  also  measured  and  the 
results  agreed  well  with  the  previous  measurements  of 
Farabee  et.  al.  LlO].  However,  the  burst  rate  computed  with 
an  Emmons  type  spot  growth  model  was  shown  to  be  inaccurate. 

When  a  small,  two  dimensional,  downstream  facing 
step  was  put  in  the  plate  near  the  leading  edge,  most  of 
the  statistics  of  the  pressure  field  did  not  change.  The 
only  exception  was  the  longitudinal  cross-spectral  density 
magnitude.  With  the  step  in  place  this  function  showed  about  a 
L  5  times  faster  decay.  The  step  also  decreased  the  trans¬ 
ition  Reynolds  number  from  about  2  x  10®  to  1.2  x  10® 

(based  on  x ) . 

Further  studies  of  the  intermittent  region  of  the 
boundary  layer  are  warranted.  The  study  of  surface  roughness 
effects  should  be  extended  to  include  several  types  of  two 
and  three  dimensional  roughnesses  commonly  found  in  marine 
and  aircraft  structures  (i.e.  seams  and  rivets). 
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An  improvement  of  Emmons'  [9]theory  to  account  fully  for 
overlapping  bursts  in  the  computation  of  the  burst  rate  and 
burst  period  would  be  valuable.  Precise  measurements  of 
the  turbulent  wall  shear  stress  in  the  intermittent  region 
would  also  be  useful.  This  would  aid  in  the  prediction  of 
the  wall  shear  stress  for  low. Reynolds  numbers. 
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APPENDIX  A 


COMPUTER  PROGRAM  LISTINGS 


A.1  DATA  A0UISITI3N  PR33RAN  !F0RTRAN  CODE) 


:ON3ITI3NALLY  A VERA3BD  CROSS-SPECTRAL  DENSITY  PR33PA* 


THIS  PR33RAN  IS  DESIGNED  TO  SUPPLE  DATA  AND  THEN 
CONDITIONALLY  REJECT  OR  KEEP  THE  DATA  DEPENDING  ON  A 
PRESET  TRIGGERING  LEVEL. 

CHANNEL  I  IS  THE  TRIGGER  III  TWO  CHARNEL  NODE,  AND 
CHANNELS  X  AND  Y  ARE  THEN  INPUT  INTO  THE  CROSS-SPECTRA! 
ANALYSIS  PROGRAM  OF  RICH  DEJONG. 

OUTPUT  INCLUOES  THE  RN3  BEFORE  AID  AFTER  CONDITIONING, 
AND  THE  CROSS-SPECTRAI.  OUTPUT  FOR  CPLOT. 


THIS  IS  VERSION  E 


NOV  CHARLES  3EDN FT 


I NPLICI T  INTE3ER*2  (I-N) 

OINSNSION  P  3X1  256)  ,»3  Y!  25  6 )  ,  CO!  256)  ,  OH  A  DC  255  )  ,  TA  B  (  256 ) 
DI1ENSION  SENS(2)  , SA I N! 2)  ,SCALFX( 2),NDATA (2,512) 
DIMENSION  RNS(2)  ,CRIN3(?)  ,  RN SR  A T! 2 )  ,  DRNS! 2 ) , R A  NO! 2) 
DIMENSION  33RN(2),S3RN(  2)  ,D JNNY! 4 ) , IDATA! 2, 2D4R ) 

DATA  SENS, 3 A  IN, 3 AND/ 0.0,  0.0, 0.0, 0.0, 0.0,1 50  00./ 

DATA  ISANP.MSANP,MSAHP/2DUR,512,514/ 

DATA  MDSV,N3ir3,NOUr/6,),255/ 

DATA  FREQ,SHO, A VE/15DDD. ,21.3,250./ 


INITIALIZE 


CALL  PR  I V 

CALL  SI NTABITAB, MBITS) 

1  1  HRI TE( NDEV, 12) 

12  FORMAT!*  CERTAIN  VARIABLES  ARE  PRESET*,/, 
1  /,•  TO  MODIFY ,  TYPE  NEW  VALUES*,/) 


TO  ACCEPT,  HIT  RETURN*  , 


SENS ( 1  )»0.  D 

100  CONTINUE 

HRITEI NDEV, 101) SENS, 3 AIN 

101  FORNATC  SENSITIVITY!  <  ,Y)  Al 
1 */1X , 4F1 0 .5 ) 

READ  !  NDE V, 1 02) IDURHYIT >,1*1,4) 

102  FORI  AT( ’4F1 0 . 2  ) 

IFCDUN1ICD.E0. 0.0)33  TO  113 
DO  103  1*1,2 

J*I  +  2 

SENS(I)«DUHNY(I) 

103  3AIN(I)*DUN1Y( J> 

33  TO  100 


GA  IN ( X, Y  )  (DR)  CSFNS-RE  1V/REF) 


11C  URITE(NDEV,111)FRE0,BAN0 

111  FORMAT! *  MAX  F.REJ  LD-OUTPUT  BANOS-HI 


(HZ) *,/,1X,3F1D.l) 


L 


■ 
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REAO(1DEV,112)OU4NY(T>,I»1,3) 

112  FDR1AT(3F1D.2) 

IF(DU14Y< 1) .E3.D.0)  3D  TO  113 
FRE3»D011Y<  1) 

BAND  (  1)*DtH"1Y(2) 

BA4D(2)«DUKSY<3) 

3D  TO  110 

113  INTVL*IFIX(D.2ED7/FREQ0.5) 

FRE3*0.2ED7/FLDAr(I1TVL) 

TPI3-0.0 

120  i!RirE(*DEV,121)rRIS 

121  FORIATC  TRIS3ERIH3  LEVEL  < YREF/4TCR0SEC) * /1X,F10.5) 
READ(YDEY,122)D;J44Y(1  ) 

12?  FORIATCFID.S) 

IF(DU11Y(1).E3.).D)  30  TO  123 
TRI3»D’J11Y<  1) 

30  TO  120 

123  CDiriios 

JSA1P*ISA1P 

1«C  «RirE(M0EV,141)5RD 

141  FDR1 A  T<  *  SIDOTHIRG  I H TE R V  AL *  , /,  1  X ,  F4»  0 ,  *  3A1PLSS*) 

READdDEV,  142)  DU  WY CD 

142  F0S1AT(F10.5) 

IF(DU*1Y<  1)  .E3.D.0)  30  TO  143 

sid-duiiycd 

IS«D-(IFIX(S1D)/2)*2*1 

S4D»FL0AT<IS40) 

3D  TO  140 

143  IS1D«IFIX(S10) 

TRI3»TRIG*TRI3*3!!D*1D.**<GAI!I<2)/10.>*1D.**(SE1S<2>/1D.)/PRE3/ 

1FRE3*2.58419E12 

KSBD«(I51D*1)/2 

JS1D*JSA4P-KS1D-1 

LSA1P«JS10-KS'*D*1 

KSA1P*JSA4P-1 

150  «RITE(ND£V,151)  RYE 

151  FDR4AT(1X,F10.2,*  AVERA  3ES* ) 

REXOCHDEY, 152)00141(1) 

152  FDR1AKF10.5) 

IF(DU41Y(1).E3.D.O)  30  TO  153 
A VE*DU11Y ( 1 ) 

3D  TO  150 

153  LIHIT«IFIX(AVE) 

HRIFE( 4DEV, 1 60) 

160  FORIATC  FOR  IlfTEPHirTAIICY  OILY,  TYPE  0*,/, 

1  *  FDR  SPECTRAL  AXALYSIS.  TYPE  V) 

READ(NDEV,151)IPFT 

161  FORIAT(II) 

HRI TE( NDEV, 1 70) 

170  FORIATC  TYPE  1  TD  10DIFY  VARIABLES* ,/, 

1  *  TYPE  0  TO  START  SA1PLIR3*) 

REA  D(  If  DEV  ,171)  ID 

171  FORIATC  ID 


ooo  ooo  ooo  ooo  ooo 
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tP( TD.NE.O>  GO  TO  103 


CLEAR  OUTPUT  ARRATS 

1000  03  1001  I-1,N0JT 
PSHI)-1.E-30 
PST<I>-1.E-33 
03(1)  -1.E-30 

1001  3UA0<I>-1.E-3? 

saitot-o.o 
00  1002  1-1,2 
RRS<I>«0.0 
DRNS(T)«0.0 

loo?  crnsid-o.o 

ICOUNT-O 

BUR3T-0.0 

IAIN  AVERAGING  LOOP 
A1K-0. 0 

1003  ICOUNT-ICOUNTH 

IF(  tC0UMT.3T.LIHT)3O  TO  1200 
00  1004  1-1,2 
S3RN ( I ) -3.0 

1004  53RKD-0.0 

AVERAGING  LOOP  FOR  DATA  COLLECTION 

110C  I3AN-0 
ITER-0 
LERR-0 

READ  DATA 

1101  CALL  CONVRT  <IDATA,3,I3ANP,INTVL,LERR) 
IF(LERR)91,1103,91 

9  1  CONTINUE 

KRITE(N0EV,1132)LEPR 

1102  F3R1ATC  A/D  FAILED*, 120) 

GO  TO  1300 

1103  CALL  UNPACK <IDATA,ISANP) 

DO  210  1-1,2 

NSUfl-0. 0 

DO  200  J»1,ISANP 
200  NSUR-NSUR*IDATA(I,J>/ISANP 
DO  210  J-1 ,IS A1P 
210  IDATA(I,J)-IDATA(I,J)-N3U1 
IFL3-1 

DEFINE  INTERRITTANCT  FUNCTION 

DO  1112  J-KSNO, JS10 
OETEK-O.O 
03  1109  I-1,ISN0 
1-J-KSRO+IO 

IDIF-IDATAC2,R*1)/4-IDATA (?,N-1)/4 
DIFfT-FLOAT(IDIF) 


UOO  UOfl  o  CM) 
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DIPFT»DIFFY*DIFFT 
DETEK-DETEK+DIFFT 
1109  CONTINUE 

IFCDETEN.LT. 1913)50  T9  1111 
IF!  IFLS.LT.UBURSr-BURSTd.O 
IFL3-1 
I3A1-I3ANO 
00  1110  1-1,2 
ADAT-FLOAT! IDATA  ( I ,  J )  ) 
S0R1<I)-S0RN(I>*ADAT*ADAT 

mo  continue 

30  TO  1112 

1111  IFL3-0 

00  1135  1-1,2 
ADAT-FL0AT<IDATA(I,J>> 
33RN(I)-S0RN<I>«>A3Ar*A0AT 
IDATA(I,J)»3.0 
1135  CONTINUE 

1112  CONTINUE 

IS  DATA  ARRAY  FILLED  t 

ICEND«ICSN0>NSA1P-1 
DO  1115  1-1,2 
D3  1115  J«ICS13,<END 
JJ-J-KSNO+1 
1115  NDATA(I,JJ)«IDATA<I,J> 

AKK-AKK*1.D 
T  TP RaTf PR49 

IF(IGA1.GE.NSANP)30  TO  1123 
IFdTER.ST. 100)30  TO  1121 
30  TO  1101 

1121  «RITE<NDEV,1122) 

1122  FORMAT! *  3ANHA  TOO  LOB*) 

30  TO  1500 

FIRD  CONTIONAL  RNS  RATIOS 

1123  3AN*FL0AT(I3AN) 

SANTOT-SANt-SANTOT 
SANP-LSAfIP-ITER 
DO  1124  1-1,2 

RNS(I)-(S0RN<I)tS0RN(I))/SANP*RN3<T) 

CR.N3(I)*S0RM(I)/3AK*CRN3(I) 

DRHS<I)»S38N(I)/<SA!lP-3MI1*1.E-33)*DRNSm 

IF((SA1P-GAS).LE.0.3)DR1S(I)-0.3 

1124  CONTINUE 

IF( IFFT. NE.O) 30  TO  1125 
30  TO  1303 

CALCULATE  FOURIER  TRANSFORM  AND  ADD  TO  RESULTS 

1125  NEI3HT>(SANP/(3AN*1.E-33))*».25 
CALL  FFTCfNDATA, TAB, MBITS, 1) 

CALL  CFS  (NDATA(1,1),X1) 

CALL  CF3(NDATA(2,1),f1) 

PSX(  1 )-PSX( 1)*4#*X1*X1*NETGHT*BEI5HT 


UIIIIO  0(10  (KKIIIII 
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psn  i)«Psr(i)*i».*Ti*r  i*aEi3m*HEi;«T 
coo  )«CD(i>*4.*r  iui-ueishmiieisrt 
DO  1126  J*2,N0UT 
:ui  CFS< MDKTKC 1  ,J),X1) 

CALL  CF3(*DATA<2, J), t 1) 

JJ*1SAHP-J 

CALL  CFSf RDA TA(1,JJ),X2) 

CALL  CFS(NDATA(2, JJ),Y2> 

XR*<XUX2>*4EI3ir 

XI«<  Y1-Y2)*HEI3-ir 

YR*<YW2)*WEI3HT 

YI*<X2-X1)*UEI3RT 

PSX(J)»PSX<J)tXR*XR*XI*XI 

psy<  J)*psy<  j)*n*YR*ii*ri 

C0<  J>*CO(J)«-YR*CR*YI*XI 
1126  30A0( J)*3UAD( J)-rR*KI*ri*XR 
30  TO  1003 


interiittancy  ootpot 

120C  NRITE(NDEV,1201) 

1  FORHATC  RMS  UAH, T3RB, TOTAL)  *> 

FACTS*. D001525/33RT(AVE) 

DO  1202  1*1,2 
RflS(I)*S0RT(RH3CI))*FACTR 
CRNS(I)*S3Rr<CR1S(I))*FACTR 
DRHS(I)*SORT(DRXS(I))*FACTR 
S138ATCI)*CR1S([)/DRHS(r) 

1202  4RITE<SDEV,1203)DRHSf I),CRR3<I),RHS(I) 

1203  FORNAT(H,3F10.5> 
SAHTOr=FLOAr(LSA!IP)*AXX 
SAH-GAHTOT/SAHTOT 
VRITE(KDEV,1204) 

1204  FORHATC  R«!3  RATIO(l,Y)  3A4RA  • ) 

HRITE(N0EV,1203)RHSRAT(1),RHSRAr(2),SAH 
BURATE*BURST*FRE0*2./SA!iror 
BURPES*3AH/(  BORATE  •fl.E-30) 

IF( BORA TE.E3.O.0)B0RPE9 *0.0 
»RITE(NDEV,1235> 

1205  FOR1ATO  AYE  BURST  RATE,  PERIOD* ) 
URITE<  If  DEV,  1206)  BORATE,  30  RPER 

1206  FORHAT(1X,2E10.3, *  (SEC)*) 

RETURN  TO  BEGINNING  TO  RESET  VARIABLES 

IF( IFFT.E3.0)  30  TO  100 

OUTPUT  FOR  SPECTRAL  ANALYSIS  SECTION 

CALCULATE  SCALING  FACTORS 

FAK*S.*5.*256. 

SCALE*4.34294*AL03CFAX/(FRE3*AVE) ) 
SCALEXC 1 >*SCALE-SENS< 1>-3 AIN (1 ) 
SCALEX(2)*SCALE-SENS(2)-3AIN(2) 

SCA LE*0. 5*( SCALSK ( 1)*SCALET(2)) 


C>  CJ  «»  IIIIU 
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PRINTED  OUTPUT  3R  PUNCHED  TAPE  t 
1300  HRITE(1DEV,1301> 

*301  FORNATC  CHOOSE  PRYNTCOS  )  ,P3NCH<  3 1 )  ,  RESTARTt  00)  , STOPC  - 1 )  •  ) 
mD(VDEV,1302>  ICH 

302  FORMAT!  1 2 > 

IF(ICH.SD.O)  33  T3  130 
IFtICH.LT.3)  33  TO  1500 
HRITE(NDEV,1333) 

303  FORMAT!  *  FRE3  PSD-X  PSD-r  CPSD-XY  PHASE  C3HER* #/) 
DF»FPE3/255. 

1L*IFIX!BAND!1  )/DF*1.5> 

NU*IFIK!BA*D!2)/DF*1.S) 

IF(R0.Sr.R3(JT>!1J«RO0r 
03  1304  1*41, *3 
FR*DF*FLOAT!I-1) 

PSK0UT*4.34294*ALDG!PSX!T))+SCALEX!1) 
PSI3UT*4.34294*AL0G!PSYf I))*SCALEX!2) 
XHA3*CO!I)*CO!I>*3UAO!I)*OUAD!I> 

C0HER*XMA3/!PSX!I)*PSY!I>> 

XHA0*2.17147*ALD3!XMA3*1.E-30>*SCAL£ 

PH  *ATAN2!0DAD(n,  CO!  I )) *57.29573 
130  4  HRITE!  ICH, 1305) FR.PSKOUr.PSYOOT.IHAG.PH.CORER 
1305  FORMAT! 1X,F6.0,4F7.1,F5.2) 

3D  TO  1300 

3Uir  PROGRAM 

1500  CONTINUE 
3T3® 

SNO 


A. 2  DATA  AJUISITIO!  SUBPROGRAMS  IINTBR3A TA  ASSEMBLY  CODE) 

•SUBROUTINE  UNPACK 

•ARGUMENTS  RDATA  DATA  ARRAY 

•  NSAMP  NO.  SAMPLES 

•UNPACK  CONVERTS  STORAGE  OF  N3ATA  FROM  BYTE  T3  HALFHORD  FORM 
2NTRT  UNPACK 

UNPACK  STM  RX .SAVE  SAVE  REGISTERS 
LH  NDATA,2(15)  FETCH  ADDRESSES 
LH  MS A  HP#  4( 15) 

LH  NSAMP#0(N3AMP)  LOAD  CONSTANT 
•SET  POINTERS 

AHR  NSANP, NSAMP  NSAMP*NSAHP* 2 

SIS  NSAMP, 1 

LHR  1 1D.NDA  TA 

AIR  MID, NSAMP 

AHR  NSAMP, NSANP 
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LHR  top, rda  ta 
AHR  T3P,NSA1P 
•JRPACK  DATA 
LOOP  19  TEMP.r»C*«ID) 
sxa*  teip,teip 
stk  tenp,o(T3P) 

313  T0P,2  BLIP  POUTERS 
313  110,1 

clhr  top,rdata 

31 L  LOOP 

LH  Rf , 3A  V  E  RESTORE  REGISTERS 
3  5(15)  RETURN 
SAVE  03  13 
R<  EOJ  11 
YDATA  E30  11 
N3A1P  E3U  12 
RID  EJ’J  13 
TOP  E3U  14 
T31P  S3U  15 
310 

ENTRY  PRIY 

*  CALL  PRIY  TO  ALLOW  I/D  INSTRUCTIONS  PRIY  REEDS  TO  BE  CULLED  ORLY 

*  DNCE  PER  PR33R3H(R3  INLINE)  EXECUTION,  RU*  ADDITIONAL  CALLS  DON'T  HURT 

*  ASSURES  FORTRAN  CALL  S3  R11  ♦  14  CAB  ««E  9 ASHFD 

PRIY  LI  14 ,X*9C*  SET  3L0  SYC  3  PARRS 
LHI  1 1 ,PRG  SET  POINTER  TO  PR 3 
STH  11, <'93' 

SVC  3,0  SWAP  PSW’S 
PRC  STH  1 4 ,X ' JC  *  RESTORE  DLO 
9  2(15)  RETURN 
END 

TITLE  UNPACK  DATA  FR3R  A3C 
ENTRI  UP3ARP 

*  UPSARP(IA,NSARP) 

*  IA  IS  INTE3ER*2  ARRAI  FR3R  CRYRT 

*  NSA1P  IS  1UR8ER  3F  3ATA  SARPLE'  t A XFN 

*  DATA  IS  UNPACKED  TO  ARRAI  FOSSA T  (2,NSA*P 
CiPSAHP  STN  9 ,RS A  2 

LH  N,  4  ( 1 5 ) 

LH  K, 3(1) 

LH  AADR, 2 ( 15 ) 

AHR  1 ,  N 
AHR  AADR  ,N 
LH R  3A08, AADR 
AHR  B ADR ,1 
LIS  TWO, 2 
SHR  B ADR , TWO 
3HR  AADR, TWO 
LHI  LI ,L31 
LIS  FOUR, 4 
L31  LB  01, 0( AADR) 

LB  D2,1( AADR) 

EX  B  R  D1,3 1 
SIBR  D?,D2 
STH  D1 ,-2 (BADR ) 
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3TH  02,0(BADR) 
3HR  AADR,T40 
SHR  3ADR, POUR 
SHR  R,T40 
3PB  LI 
L.1  R,RSA»E 
AH  15,0(15) 

33  15 

R3AFF  05  16 
1  E3U  B 
AADR  E30  5 
8ADR  E3U  10 
T40  I3U  11 
FOUR  EJU  12 
01  EOJ  IB 
02  E  Q'J  14 
LI  EQJ  15 
3RD 


a  .3  subprogram  *33Rm*  ( irterdata  assembly  oode) 

TITLE  ****•**•*••**•  SAMPLE  a/D  CORYFRTER  •**•**•*•**• 
ertri  comr  (a,ichai,rsaip, period, lerr) 

*  FORTRAR  CaLLaBLE  3J3R00TIRE 

*  a  13  INTEGER* 2  A( 2,1SAR»/2)  ARRA  OF  DATA  (RCHAR*0) 

*  IRT3SR*2  A(R3ABP/2)  IF  RGHAR  ROT  ZERO 

*  PERIOD  IS  IRTF3ER  3IVIR3  RUBBER  OF  HIGROSECOROS  6ETHEER  CORYERSTORS 

*  RSA1P  IS  RUBBER  OF  rORYERSI ORS  TO  PE  MADE 

*  LERR  IS  L03ICRL  V RRIR 9LE 

*  3ETURRED  TRUE  IF  ER03  033URS,  ELSE  .FALSE. 

*  RCHRR  IS  2HARREL  TO  SAMPLE 

*  D  BOTH  CHARNELS.  FIRST  DATA  4 ILL  BE  3RA11EL  A 

*  1  CHANNEL  a  ORLT 

*  2  CAHRNEL  9  ORLT 

*  EYE3Y  CONVERSION  TRIES  OP  ORE  BITE  OF  SPACE 

*  output  oara  is  biploar,  (b  3it  sisrbo  scaled  fraction) 

CDNVRT  sY«  0,RSAV£ 

LHI  11,1*20*  SET  UP  331 HARD 
LIS  10,1  T40  BYTES  PER  SAMPLE 
LB  12,2(15)  3ET  ADDRESS  OF  A35S 
LHR  3,12  SET  BS3IRRIR3  OD  DATA  AREA 
LH  13,0(13)  SET  RCHA1 
BZS  BOTH  DOIRS  BOTH  CRARREL5 
OHR  11,13  OR  IR  CHARNEL  ROBBER 
OHI  11,2  TURR  OR  ORE  CHANNEL  BODE 
LIS  10,0  ORE  BITE  PER  SAMPLE 
BOTH  BOO  • 

LB  14,0(14)  SET  RSAHP 

SLHl  14,0(10)  MULTIPLY  Bf  LERSTH  OF  ITER 
AHR  12,14 


THIS  ?kOt  is 
TROW  ODi  t  * 


„K,  pun  n*a*»* 

:  T0  ^  


(JU(M)IIIMI  DO 
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AIS  12,1  GST  END  OF  S'JfFER 
LHR  4,12  SET  END  OP  OR TA  R REA 
LH I  1,X*S6* 

OC  1.ST0P  STOP  IT  FIRST 
LH  12,0(15) 

313  12,1 
HHR  1,12 
LI  3  12,0 

SPSB  0,12  FORCE  0  STATUS,  S*FE  OLD  STATUS 
OCR  1,11  OUTPUT  START  CONN  AND 
RBR  1,3  READ  THE  DATA 

EP3R  12,0  RESTORE  OLO  STATUS,  SAVE  ERROR  FLAG 
OC  11, STOP  STOP  IT  AGAIN 
LHR  12.12  TEST  ERROR 
32  S  OK 

LC3  12,1  -1  •  FORTRAN  .TRUE.,  0». FALSE. 

0<  LH  1S.RSAVEOO 
LH  15,10(15) 

3TH  12,0(15)  STORE  ERROR  P»  RANTER 
LN  O.RSAVE  RETURN  TO  FORTRAN  PROGRAN 
AH  15,0(15) 

3R  15 

STOP  DB  X*10* 

03  * 

R3AV  E  D3  32 
2ND 


A  .4  PRESSURE  STATISTICS  PLOTTING  PROGRAN  (FORTRAN  CODF) 


PRCGRAN  TO  PLOT  THE  HALL  PRESSURE  RTA  TISTTCS  NEASURED  IN  A  BOUNDARY  LAIPR 

THIS  PROGRAN  INPUTS  POUER  AND  CROSS  SPECTRAL  DENSITIES  ?P0" 

THE  DATA  AJOISITION  PROGRAN,  N0N-DINENSI0NALITE3  THEN  AND 
PLOTS  THEN. 


INTEGERS  XL, 10 


OINENSION  ARRAK  2,5,201)  ,XS(4)  ,KL  (  4  0  )  ,  ID(  8  )  ,  A  (  1 1 , 201  )  ,  «»ET(  2  ) 
1  ,R0'JGH(2),NVAR(  11) 

LOGICAL  KNAR 

DATA  ROUGH  /3 H3N0 , 3HR 00 / 

DEFINE  FILE  1 D( 20, 1250, J ,NRP) 

DO  5  I«1,11 
5  NTAR(I)«. FALSE. 

INPUT  DATA 

GENERAL 


REA0C3, 10)NRUN5 
10  FORNAT( 13 ) 

DO  190  IT»1,NRDNS 
TAPE  DATA 
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READ<B,20HID(L)  ,  L*1,8)  , NRP , L0C,N PLOT ,NPR ,NN A* , LR 
20  F0R1AT(9A2/6I7> 

C  TEST  CONDITIONS 

READ(8,3D)US,0L,RS,RL,D?,DPH,TS,3ANS 
30  F3RNAT<9F10.3> 

C  FIND  DATA 

S3  TO  (40,130), LOC 
C  DATA  ON  FILE  (ASCII) 

40  DO  60  K*1,2 

C  READ!  12, 53)  HARR  AT  (K,  I,  J>  ,1*1 ,6)  ,  J-1,  NNAY  ) 

READ(B,55>((ARRAY<K,I.T),T*1,6),J*1,NNAY) 

C  50  F3RSAT(1X,F7.0,:iF7.1,F6.2> 

55  FORHAT(1X,F7.3,4F7.1,F5.2) 

50  SON TI IDE 
c  :a  LCD  la  te  oensitt 

R0*-0.0023*TSM.352 
C  STATISTICAL  calolatidns 
C  INITIALIZE 
HPH*  1 13 
8RCH»D.D 
LFLA3*0 

3S*R0»US*US*0.5 
RSDJ«RS/OS 
RLDU»RL/UL 
RS0UP*RS0U*57. 29578 
DSDU*DS/US 

PSDIN*0.0795798/D300/OS/QS/S3RT(3A85) 

C  SPECTRAL  CALCULATIONS 

IF(  ARRAK  1,5,1 >  .  LE. D.  0)  9RCH*BRCR«-36  D  .  D 

IF(  ARRAK 1,5,1>-ARRAY(1,5,2>.3F.DPH)LFLA3*1 

DO  90  J*1,<mY 

J1*J+1 

J2«J*2 

N*ARRAY(1,1,J)»S. 28319 

C  ADD  2  PI  TO  PHASE  FDR  EACH  BRANCH  OF  ARCTAN  CROSSED 
ARRAK  1,5,  J)*ARRAY(  1,5,  JKBRCH 
C  ADD  ANOTHER  2  PI  IF  CROSSING  A  BRANCH 
BRCH-3RCH*35D.D*FL0AT(LFLA3) 

LFLASO 

IF( J.3E.NNAX-2>30  TO  7D 
C  TEST  FOR  BRANCH 

IF(ARRAY(1,5,J1)-AFRAY(1 ,5, J2).3E.DPH  )LFLAG*1 
I F< ARRAY <1,5, J2)-ARRAT( 1 ,5, J1 ) .3E.DPH  )LFLAS*-1 
C  CALCULATE  NDN-DIHENSIDNAL  STATISTICS  AND  ERRORS 
70  A(1, J)«H*RSOUP/ARRAT< 1,5.J) 

A(2,J)=H*DS0U 

A( 3, J)*N*RSOU/A( 1 ,  J) 

A(4,J)*N*RL0U/A(1»J) 

DO  80  It*  1  #2 
N1-K+4 
N  2*K  +6 

ARRAKK,2,J)»0.5MARRAY(X,2,J)+ARRAKir,3,J>> 

A  (  N1  ,J)a2»D*(ARRAY(K,3,J)-A9RAY(!C,?,J)) 
A(.N2,J)a10.**((ARRAY(IC,»,  J> -ARRAY (K,2 , J)  3/1D.0) 

B  0  3311  tlliOC 

A(9,J)a1D.0*ALD310<1D.*»(ARRAY(1,2,J>/10.O>*PSDIR> 
A(1D,J)*ARRAY(1,1 ,J) 
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11(11,  J)-ARRAT(1,5,J> 

9 c  continue 

UR  ITE  N3N-DINENSI3NAL  SPECTRA  31  FLOPPY  DISK  (Bid) 

URITEC  13*NRP)((A(I,J),T«1,11),J=1  ,MAX> 

33  TO  110 

READ  K3NDINSN5I0NAL  SPECTRA  FR3H  FLOPPY  DISK  (BIN) 

130  READ(13*NRP)( (A(I,J),I*1 ,11),J*1,4NAX) 

PR  Id  T  NJHBERS  DN  LP  T 
110  30  TO  < 120, 150) ,NPR 

12  0  URITE(5, 130)  (  T3(L),L-1,R  )  ,0 3  ,  RS,  TS , RO , DS,  OL,  P L,  R 3U3H (  L*> ) 

13C  FORK AT( 1 H 1 , 0 A ?/1 2 X , * U IN •  SEP . * ,T33, * TFNP  3F1S  DEIST 

IV  (STR)*, 4X, 5513. 3/*  (LAD*  ,4K,  2E1 0 ,3/1  OX  ,  A  3  ) 
URITE(5,143)((A(I,J),I=1,11),J=1,4NAX) 

140  FORNATCO  'JC/J  W  .D5/U*  ,  T 29  ,  •  W  .R /fJC  *,  T4  3  ,  *  ERROR  t>HT<R)/PH 

11(H)  PHinn  FRE3  PHASE*/*  ( STR  )  (STR)  ( 3TR ) 

2( LA  T )  (STR)  (LAT)  (STR)  <L»T)  (43d -DIR)  (HZ)  (DEO ) * / 
3(F3.4,3F9.3,F5.1,F7.1,2F9.4,F9.1 ,2F13.1)) 

HAKE  PLOTS 

1 5 C  IF(NPL3T.E3,D)33  TO  190 
D3  155  J="PH,N1AX 
A(1,J)«3. 

A ( 3, J ) »D . 

A(4,J)*3. 

155  CONTINUE 

DO  1 B 3  IP  =  1 ,NPL3T 
READ(3,150)NV,I3,NX,XS,(L 
160  F3RNAT(3I7,4F7.3/40A2) 

03  170  L= 1 , 3 
Xl(L)*I3(L) 

1 7 C  CONTINUE 

NVAR(NV)».TRU5. 

CALL  PICTR(A,  11,XL,XS,NXAR,»!!AT,NX,-l,  1334  ,rs,FT,1 ) 

NYAR(NV)*. FALSE. 

PA'JSE 

19C  CONTINU6" 

190  CONTINUE 
CALL  EXIT 
END 
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APPENDIX  B 

CALCULATION  OF  THE  THEORETICAL  INTERMITTENCY 
AND  BURST  RATE  FUNCTIONS 


Emmons  [9]  derived  the  following  relations  for  the 
intermittency  and  burst  rate  at  point  P(x,y)  on  a  flat  plate. 


Y  (P) 

fB(p) 


=  i-exp[-  g(PQ)dV  ] 

P 


r  g(pQ)  r 

=  -  exp[-  g(P')  dV']dV 

JR  At  JR'  °  o  o 


(B.l) 


(B.  2 ) 


where  the  point  P0  is  in  the  "retrograde  cone"  R  and  the 
point  Pq  is  in  the  "truncated  cone"  R'  in  the  (x,y,t) 
coordinate  system  (see  Figure  B.l).  g(PQ)  is  the  source 
rate  probability  density  function  and  At  is  the  duration  of 
the  burst  at  point  P. 

The  integrals  in  equation  B.l  and  B.2  can  be 
solved  if  the  volume  element  dVQ  is  written  (assuming  a 
triangular  burst  shape) , 

dvo =  r-  <*-*o)2av  (B-3> 

00 

u  Ua> 

where,  o*  =  °°  tana 


(B.  4) 


FIGURE  B.l  Rectroqrade  and  Truncated  Cones  Shown  in  the  x.y.t.  Coordinate 
System  After  Emmons  [9].  (A  Triangular  Turbulent  Spot  Shape 
is  Assumed). 


frees tream  velocity 


half  angle  of  spread  of  bursts 


dV  =  A tdx  dy 


Case  Is  Line  Source 


The  first  assumed  form  of  g(P  )  is  Dirac's  Delta 


function 


where  n  is  the  number  of  sources  per  unit  length  per  unit 
time  along  the  line  x  =  x  .  Substitution  into  equation  B.l 


yields 
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Where  the  Heaviside  function,  H(x)  is  used  since  the  value 
of  this  integral  strongly  depends  on  the  relative  locations 
XQ  and  x^.  Equation  B.2  becomes 


ny+d  t 

y-d  "S(,<0  '  Xt)eXpt’  it  <X  '  Xt>2,<<Xo  *  Xt)JdyodX< 


(B.  9 ) 


where  d  =  (x  -  x^)  tana  (B.10) 

Solving  the  integrals  in  equation  B.9  gives 

fB  =  2ntana(x-xt)exp[-  ^  (x-x  )2]  .  (B.ll) 

00 

The  nondimensionalization  of  ffi  is  done  as  in  Farabee  et.  al 

[10]: 

fB*  *  V" 

Combining  equations  B.7,  B.ll  and  B.12  yields, 

f B*  =  2  /(l-Y ) In  (j—)  (4.12) 

and 

V  ■  °'420  fB  4V°.  •  (4.13) 

Case  2 :  Constant  Source 

The  second  form  of  g(P  )  is  a  constant  for  x  >  x  . 

<-»  o  —  t 

°  <XQ  <  V 

g(Po>  { g  (X  >  X  )  (B.  13) 

o  —  t 


tana 


(B.12) 


-110- 


Equation  B.l  now  becomes. 


(P)  =  1  -  exp[-  (x  -  xt)3] 


(B. 14) 


and  B.2  becomes 


fB  =  2  tan  a  g  AXT2  f  nexp[|^3-  AX.J,3  (n  3-x3)  ]dn 

J  O  oo 


(B. 15) 


x  -  x. 


where  x  = 


(B.16) 


axt  =  x|y=0>99  -  xlY=o.oi 


(B. 17) 


Combining  B. 14  and  B.17  yields. 


AX„  =1.45  3 

A  il 


(B. 18) 


Rewriting  equation  B.15, 


U.  2/3  fx  3 

=  8.72  gtana  ( - )  nexp[3.04 (n  -x  )]dn. 

go*  Jo 

Nondimensionalizing  fn  again  as  in  Farabee  et.  al.  [lO], 


(B. 19) 


fB*  5  fB/gtana  ( - ) 


(B. 20) 


Combining  equations  B.18,  B.19  and  B.20  gives  the  desired 


relations, 


fB*  =  8.72  jXnexp[3.04(n3-x3)]dn 


(4.15) 


